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Abstract 

We present an exact solution to the problem of the relativistic motion of 2 point masses in (1 + 1) 
dimensional dilaton gravity. The motion of the bodies is governed entirely by their mutual gravitational 
influence, and the spacetime metric is likewise fully determined by their stress-energy. A Newtonian limit 
exists, and there is a static gravitational potential. Our solution gives the exact Hamiltonian to infinite 
order in the gravitational coupling constant. 
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The problem of motion is a notoriously difRcult one in gravitational theory. Although approximation 
techniques exist jl], in general there is no exact solution to the problem of the motion of N bodies each 
interacting under their mutual gravitational influence, except in the case N = 2 for Newtonian gravity, or in 
(2 + 1) dimensions, where the absence of a static gravitational potential allows one to generalize the static 
2-body metric to that of two bodies moving with any speed . 

We present here an exact solution to problem of the relativistic motion of 2 point masses under gravity 
in (1 + 1) dimensions. The dimensionality necessitates that the gravitational theory we choose is a dilaton 
theory of gravity; however our choice of dilatonic gravity is such that the dilaton decouples from the classical 
equations of motion ^ . Consequently the motion of the bodies is governed entirely by their mutual gravita- 
tional influence, and the spacetime metric is likewise fully determined by their stress-energy [Q. Unlike the 
(2 + 1) dimensional case, a Newtonian limit exists, and there is a static gravitational potential. Our solution 
gives the exact Hamiltonian to inflnite order in the gravitational coupling constant. We can thus view the 
whole structure of the theory from the weak fleld to the strong flcld limits. 

We begin with the A/'-body gravitational action 
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where 4* is the dilaton field and g^^, g, are the metric and its determinant, R is the Ricci scalar and is 
the proper time of a-th particle, with k SvrG/c'*. It is straightforward to show that the system of field 
equations reduces to 

and 

iv^W.^- - Qv^Wa* - V^'J'^ - V^V,* = kT^, (3) 
where the stress-energy due to the point masses is 

T^,u ^^^^a j dTa-^g^ag,.p^^S^{x - Za(Ta)) , (4) 

and is conserved. Note that is a closed system of + 1 equations for which one can solve for the single 
metric degree of freedom and the A^ degrees of freedom of the point masses. The evolution of the dilaton 
field is governed by the evolution of the point-masses via (|^). It is easy to show that the left-hand side of this 
equation is divergenceless (consistent with the conservation of T^^), yielding only one independent equation 
to determine the single degree of freedom of the dilaton. 

We shall work in the canonical formalism for which the action (|l|) is written in the formg 

J dx^i^Y.PaZaSix-Za{x"))+7rj + Ilii + NoR'^ + NiR^^ (5) 
where 7 — gn, No = {—g^^)~^^'^ , A^i = gia, tt and 11 are conjugate momenta to 7 and ^E" respectively, and 



(6) 

i?i = ^TT- in*' + 2^' + y . (7) 
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with the symbols ( ' ) and ( ' ) denoting and di , respectively. 
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The quantities Nq and A^i are Lagrange multipliers which enforce the constraints R'^ = = R^. Since 
the only linear terms in these constraints are (^''/^/7)' ^i^d tt', we may solve for these quantities in terms 
of the dynamical and gauge {i.e. co-ordinate) degrees of freedom. These latter degrees of freedom may be 
identified by writing the generator (which arises from the variation of the action at the boundaries) in terms 
of the former quantities, and then finding which quantities serve to fix the frame of the physical space-time 
coordinates in a manner analogous to the (3 + l)-dimensional case 0, 

Carrying out this procedure, we find that we can consistently choose the coordinate conditions 7 = 1 and 
n = 0. Eliminating the constraints, the action (^ then reduces to 



(8) 



where the reduced Hamiltonian for the system of particles is H = J dxTL = J dxA^ , where A = d^/da 
Here \I' = 'i'{x, Za,Pa) and is determined by solving the constraint equations which are now 



(*')^ + l^'^T^^ + K 2_. \^Pa + ml6{x - Za) ^0 
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An exact solution to these equations in the 2-body case may now be obtained as follows. Consider first 
the case Z2 < zi, for which we may divide spacetime into three regions: zi < x {{+) region), Z2 < x < zi 
((0) region) and x < Z2 ((-) region). Writing 5* = — 41og|0| and tt = x' we find in each region that (p is the 
sum of growing and decaying exponentials in x and that x is a sum of terms linear in x. Matching these 
solutions at the boundaries x = zi and a; = Z2 of each region allows a determination of the coefficients of the 
exponentials in the -I- and — regions in terms of those in the region. 

Because the magnitudes of both and x increase with increasing |a;|, it is necessary to impose a boundary 



) vanish. This condition 
Incorporating this into 



condition which guarantees that the surface terms obtained in passing from (^ to 
has been shown to be 4"^ — 4k^x^ = 0, which must hold in the + and — regions 
the matching conditions allows one to fully determine the coefficients of the exponentials in the region in 
terms of the momenta and positions of the point masses. 

Solving the constraints and the equations for A^o and A^i yields 
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M2 
A{Mi 



M2 



PO 



Ml 



1/2 



L + — 21) 



2 Lo 
f{zi + Z2)LiL2)/Lo, e2 
1/2 



1, and 



L2_ 

M2 



U/2 



^0 



L2 

Ml 



1/2 



ill e%Lo{x-z,) ^ ±n_ 



Li 
M2 



1/2 



L_(a;-Z2) 



^-fLo{x-Z2) 



as the solutions for the relevant field variables. For convenience we have defined the quantities 



Li=Lq- M2 



L±=AX± e{pi 
L2 = Lo- Ml 
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where —7711 = r]22 = ^,'^12 = V21 = and Lq — AX — e(pi ~P2)- These quantities are related by the equation 



L1L2 = AfiMae^^o^^i-^^) (18) 

which determines X. 

The parameter e is a constant of integration associated with the metric degree of freedom. We have two 
types of solutions corresponding to e = ±1. Under time reversal, these solutions transform into each other, 
ensuring invariance of the whole theory under this symmetry. It is straightforward to show from this solution 
that the Ricci scalar vanishes everywhere except at the locations (21 (t), 2:2 (i)) of the point masses. 

A calcuation shows that the Hamiltonian of the system is H = AX. Consequently ( |l8|) gives an exact 
solution for the Hamiltonian H in terms of the co-ordinate and metric degrees of freedom. It is straightforward 
to show that the total momenta pi + p2 is conserved, and so we may choose a frame of reference with 
Pi — —p2 = P- Hamilton's equations then yield 

p = -^ALiL2 (19) 

= -^+^-^ (21) 

as the dynamical equations for the 2-body system coupled to gravity. 

Repeating this analysis for zi < Z2 yields the general equation for the Hamiltonian 



iyPi + ml- ep2 - H){^pI + ml + epi - H) 

= (v/??^^-ePi)(VPi^^ + eP2)e^^^-'(Pi-?^^»l''l (22) 

where r = zi — Z2 and Pa = PaSgn(zi — Z2). 

Equation ( |2^ ) describes the surface in (r,p, H) space of all allowed phase-space trajectories. Since if is a 
constant of the motion, (a fact easily verified by differentiation of (|2^ ) with respect to t) a given trajectory 
in the (r, p) plane is uniquely determined by setting H = Hq in (p2[). For Hq sufhciently small, the trajectory 
can be regarded as a relativistic perturbation of the Newtonian case. However once Hq is sufficiently large 
there exist a qualitatively new set of trajectories which cannot be understood as relativistically 'correcting' 
the Newtonian ones. The transition from the first set to the second set is smooth. 

In the equal-mass case, the Hamiltonian is 



H 



W (^—^ (\r\ \Jp'^ + rn? — eprj exp | {^\ \J p^ + — eprj ^ 



+ \/p'^ +m'^ + epsgn(r) (23) 



K\r\ 

where W{x) is the Lambert W- function defined via 

yey = x =^ y = W{x) (24) 

and has two real branches ^ which join smoothly onto each other. The principal branch (for which W{x) > 
— 1) in ( p3|) reduces to the Newtonian limit for small k; the second set of trajectories mentioned above are 
described by the other branch. 

Some characteristic phase-space plots for the equal mass case are given in figs. 1-2. In each of these we 
have included the corresponding trajectory from the Newtonian theory for comparison. We see that as Hq 
increases, the trajectory becomes more 'S'-shaped, with the particles reaching their maximum separation for 
some positive value of p, where the velocity f = 0. This occurs because p depends upon the momenta of 
the particles and of the metric. Under time-reversal, the trajectory for a given value of Hq is obtained by 
reflection in the r = axis. Fig. 3 shows three phase space plots for the unequal mass case, with the mass 
of particle 2 being 1/10, equal to and twice the mass of particle 1. 
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An interesting limiting case of (22) may be obtained by setting zi — z, Z2 = 0, mi = /i, m2 = m,pi = p 
and p2 = 0. Taking the lowest order terms of /i, we obtain 



H = m + yjp^ + /i^eTI^I - £p sgn(2;) |^e"™'^' - ij . (25) 

which is just the Hamiltonian of a test particle of mass /i under the influence of a static source at the origin 
of mass m, in which the metric is A^o = e'^^'^l, N\ = e sgn(2) (e^^l^l — l) [^. 

We hope to extend our analysis of this problem to the many-body problem and to more general theories 
of dilaton gravity. 
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Figure Captions 



Fig. 1 A comparison of rclativistic (s-shaped) and non-relativistic (oval) phase-space trajectories for 

H ~ 2.2m. All axes are in units of m. 

Fig. 2 A comparison of relativistic (s-shaped) and non-relativistic (oval) phase-space trajectories for 

H — 4m. All axes arc in units of m. 

Fig. 3 Phase-space trajectories for H = Ami, with m2 = 2mi (innermost curve), m2 = mi and 
m2 = .Imi (outermost curve). All axes are in units of mi. 
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